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Abstract 

From a phenomenological point of view, we study active-active and active-sterile 
flavour-changing (and flavour-conserving) oscillations of Dirac-Majorana neutrinos 
both in vacuum and in matter. The general expressions for the transition proba- 
bilities in vacuum are reported. We then investigate some interesting consequences 
following from particular simple forms of the neutrino mass matrices, and for the 
envisaged scenarios we discuss in detail neutrino propagation in matter. Special em- 
phasis is given to the problem of occurrence of resonant enhancement of active-active 
and active-sterile neutrino oscillations in a medium. The peculiar novel features re- 
lated to the Dirac-Majorana nature of neutrinos are particularly pointed out. 
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1 Introduction 



Today we have several indications in favour of non zero neutrino masses and mixing. 
The solar neutrino problem, i.e. the observed deficit of solar neutrino fluxes |]T|, is a well 
estabhshed tool whose resolution requires (almost without doubt 0) neutrino physics 
beyond the (minimal) Standard Model. The acceptable solutions to this problem, in 
terms of vacuum or matter flavour oscillations or spin and flavour oscillations P as 
well as in terms of active-sterile neutrino conversion P], all need non vanishing neutrino 
masses and mixing |0, |^, ^, |10[ . 

The second indication in favour of neutrino oscillations come from the observed deficit of 



atmospheric muon neutrinos with respect to electron neutrinos |Tl| that can be explained 
in terms of z/^ Vj. or z/^ — > z/g or even active-sterile neutrino transitions [0 . 
Laboratory direct searches for massive neutrinos only give, at present, upper limits on 
neutrino masses and the same is valid for reactor and accelerator neutrino oscillations 



experiments [|14[, except for the LSND experiment |15| whose results seem to be explained 



in terms of — z/g oscillations. 



Hints for massive neutrinos also come from cosmology, looking at Vj. as the most probable 
candidate for the hot component of the dark matter [|I6| and from the observed peculiar 



velocities of pulsars [0 . On the other hand, in Grand Unified Theories, which attempt to 
give a unified view of electroweak and strong interactions, massive neutrinos are predicted 
||T8| together with other phenomena violating both lepton numbers and baryon number 
(such as, for example, proton decay). 

However, the most intriguing fact is that a simple scenario with only three massive neu- 
trinos cannot account for the solar neutrino problem, the atmospheric neutrino anomaly 
and the LSND result. This is because the three squared masses differences Am? for the 
three oscillation solutions to these problems are all distinct between them: the resonant 
MSW solution to the solar neutrino problem requires Am? ~ 10~^ eV^, while for the 
atmospheric anomaly Am^ ~ 10~^ eV^ is needed and Am^ ~ 1 eV^ for the LSND result. 
Many analyses have been conducted for giving a unified view of the three problems 
in terms of neutrino oscillation (taking into account also the limits from laboratory ex- 
periments) and a coherent picture seems to emerge with four massive neutrinos, namely 
the three known neutrinos plus a sterile neutrino. Note that four neutrino mass eigen- 
states are needed, but not necessarily more than three neutrino flavour eigenstates. This 
scenario is easily realized if one considers neutrinos as Dirac-Majorana particles described 
by the following general mass term in the electroweak lagrangian [p!0| : 



- CT = Y.^VR M^, + \Y.^vr Mil ^iL + \ Y.nL Mf,i + h.c. (1) 

l,V 1,1' 1,1' 

Here /, /' = e, fi,T label the three flavour eigenstates and Mp, Mi, M2 are the Dirac and 
the two Majorana mass matrices which, in general, are hermitian and non diagonal (how- 
ever. Ml and M2 have to be symmetric). To construct the mass term in (|l|) we need the 
three known left-handed neutrinos (and their antiparticles) and other three right-handed 
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sterile neutrinos (and their antiparticles) After the diagonahzation of (|T]) we can obtain 
in general six mass eigenstates which are Majorana fields; so in this framework we can 
easily endow the above scenario with four massive neutrinos coming from the experiments. 
Note that if neutrinos are really described by the mass term in (|lD, the total lepton num- 
ber is no longer conserved and peculiar phenomena, as neutrinoless double beta decay 
and neutrino-antineutrino oscillations can take place. 

We stress that (|l|) is predicted in many GUTs in which the popular "seesaw" mecha- 
nism 1 20] can give rise to very small neutrino masses in a very natural way by supposing 
Ml !^ and Mo -C M2. However, this is not the only possibility; recently some models 
assuming Mi ~ M2 have been proposed ||21|] for accounting the three experimental indica- 
tions on neutrino oscillations discussed above. Here we further explore this last scenario 
and study flavour transitions of Dirac-Majorana neutrinos from a completely phenomeno- 
logical point of view, adopting no particular model. This work is a generalization to 
flavour transitions of previous papers in which we studied peculiar oscillations 

of Dirac-Majorana neutrinos. We now assume, for simplicity, only two flavours, so Mjj, 
Ml, M2 in (0) are 2x2 matrices in the flavour space. In the following section, the basic 
vacuum oscillations allowed by (|l]) are studied and transition probabilities are explicitly 
given in the general case. Some very interesting consequences due to particularly simple 
forms of the mass matrix are also investigated. In section 3, given the effective hamilto- 
nian of Dirac-Majorana neutrinos interacting with a medium, resonant matter oscillations 
are considered along with a qualitative discussion of the phenomenon with the aid of the 
level crossing diagram. Finally, in section 4, there are our conclusions and remarks. 



2 Dinamical evolution of Dirac-Majorana neutrinos 
in vacuum 



Let us consider the propagation in vacuum of Dirac-Majorana neutrinos with 4-momentum 
= {uj, k) described by the following lagrangian: 
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(for convenience we have absorbed a factor 1/2 in and Mm with respect to the mass 
terms in (|l])). In the chiral Weyl basis for the Dirac gamma matrices, denoting 
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^Obviously, the generalization to more than three famihes is possible and straightforward 
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In a more compact form, indicating 




and 



/ Md Mm \ 
= Mm M J 





the equation (|^) can be written as 

Eqs. (1^) shows that, in general, chirahty-changing transitions are possible (in fact, in 
vacuum, chirality is not in general conserved), but these are suppressed with respect to 
the chirality-preserving ones in the ultrarelativistic limit, because in this limit chirality 
almost coincides with helicity, which is strictly conserved |2^. Due to this suppression. 



chirality-changing transitions are not of very practical interest, and we will not consider 
further these processes. At the leading order in the ultrarelativistic limit {u — k 2k, 
from d^) we deduce that 

( uj - k\ M 

while the approximate equation for hl is 

HriL = ujriL (9) 

with 

H k + — (10) 

and the matrix M is given in (^). Here, for simplicity, we limit ourselves to only two 
flavours, for example e and fi, so that the mass matrices in (^ have in general (assuming 
CP conserved) the following non diagonal form: 



Md = [ ] Mm = [ "j^ ] (H) 

The 4x4 hamiltonian in (|10]) can be easily block-diagonalized by using the unitary matrix 

(12) 




where / is the 2x2 identity matrix. In fact, introducing the Majorana states rii = V n^, 
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we can rewrite (0) as 



with 



where 



H = k + 



M2 



Ml 





k + 



Ml 
~2k 



(14) 
(15) 



M± = Md ± Mm (16) 
It is now a simple task to completely diagonalize ([151) by means of the mixing matrices 



^4 



COS 0± sin 6± 
— sin 6± cos 6± 



(17) 



with the mixing angles given by 



tan 264 



2 m 
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mi 
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(m^ = m£ ± and so on). Indicating with 

= \ {^te + + V '^'^^e ~ "^^A*)^ + 4(m± )2) 



± p 



+ 4(m± )2 



the eigenvalues of the mass matrices in (|T6|), the four energy eigenvalues are 



Ef^ = k + ^ ^'^^ 



2k 



to which correspond the eigenstates 
= U 



U+ 
t/_ 



c 



where we have used the shorthand notation 



V2± 



(19) 
(20) 

(21) 

(22) 

(23) 



Given the relation ( p2D between the energy eigenstates and the flavour ones, we easily get 
the time evolution of the states created by weak interactions: 



/ I ^eL{t) > \ 

V I ^^lW > J 



V U' AUV 



( I i^eL(O) > \ 

I '^/.l(O) > 

I KM > 
V I ^^l(0) > J 



(24) 



where 

A = diag{e-^^'+\e-'^'+\e-^^'-\e-''^'-'} (25) 

The transition probabihties P{vi-^ Vj) are then given by the squared modulus of the 
corresponding matrix elements | < ^'j(O) | z/j(0) > ^ in (|2^) ; after some calculations we 



obtain 



^/ N 1/9^9 Ami 9 ^ 9 Am^ 

^'('^eL ^ i^uL ) = - sm^ 29+ sm^ — -^t + sm^ 26. sm^ — -^t 
A \ Ak Ak 

Ami Am^ S \ 
+ 2 sin 29. sin 2^_ sin — r^t sin —r^t cos — t (26) 

Ak Ak Ak J ^ ^ 

^/ X 1 / 9 ^ 9 Ami 9 ^ 9 Am^ 

P z/eL ^ = - sin2 2^+ sin^ — -^t + sin^ 2^_ sin^ — -^t 

^ A\ Ak Ak 

Ami Am^ S \ 
- 2 sin 26^+ sin 29_ sin — -v^t sin —r^t cos — t (27) 

4A; Ak Ak ^ ' 



9 9 S — Ami + Am^ ^ „ 9 S — Ami — Am^ 
Pi'^eL - = cld sin^ + clsl sm' 

,0 o S + Ami + Am^ , , , S + Ami - Am^ 
+ 4^^_si„^— .t + s^i^in^ 

- cisi sin^ _ ^2 ^2 gj^2 Am_ _^ 

and for the survival probability 

2 2 2 ^ ~ Ami + Am^ 22 2 ^ ~ Ami — Ami 

P{l^eL ^ Z^eL ) = 1 - C+C_ Sm ^ 1 - C+S_ Sm ^ 1 

2 2 2 S + Ami + Ami 22 2 ^ + Ami - Ami 

— s.c sin 1 — s.s sin 1 

+ ~ 8A; + - 8k 

2 2-2 Am^ 22-2 Am_ 

— c,s, sm — -^t — c_s_sm — -- — t 29 

^ ~^ Ak Ak 

(we have adopted the shorthand notation c± = cos 9±, s± = sin 9±). These probabilities 
in general depend on 5 parameters of the underlying theory, that is 2 mixing angles 9±, 
and 3 mass parameters Am^ = m2_|_ — m,l+, S = mfj^ + m^^ — — m2_. 
Note that in the limit of zero mixing angles all transition probabilities vanish except 
P{i^eL ^cl) which we recover the Pontecorvo oscillation formula |2^, |2^: in this 
limit only neutrino-antineutrino oscillations with no flavour change are possible with both 
Dirac and Majorana mass terms. 

The obtained results for the transition probabilities in the general case, eqs. (p!7|)-(p!9|), 
are rather complicate and then difficult to analyze. In the following we discuss some 
very interesting particular cases obtained for peculiar forms of Dirac and Majorana mass 
matrices. 
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2.1 Pure Dirac and pure Majorana neutrinos. 



For the two limiting cases 



(30) 
(31) 

(32) 



(33) 

(34) 
(35) 

As we know, in these frameworks we have no neutrino-antineutrino oscillations but only 
flavour transitions, for which we recover the standard results [0]. Note that for pure Dirac 
and pure Majorana neutrinos the result is the same, but this holds only in the ultrarela- 
tivistic limit 1^. 
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2.2 The cases of Dirac mixing and Majorana masses and vice- 
versa. 

In general, all the elements of the Md and Mm mass matrices in (0) are non zero; however, 
we can constrain these by making some physical ansatz. 

Looking at the quark sector of the Standard Model we note that quarks are Dirac par- 
ticles (obviously!) and flavour (weak interacting) eigenstates are mixed to give the mass 
eigenstates. In analogy, we can assume that also for neutrinos the Dirac mass matrix is re- 
sponsible for the mixing of the flavour eigenstates. Nevertheless, we know that neutrinos, 
if massive, are much more light than the corresponding quarks (see the limits reported in 
P^ ) so that one can think that their masses (but not necessarily mixings) are generated 
by a Majorana mass term. The most simple forms for M^ and Mm translating these two 
ansatz are then 



We now explore the implications of (pq) . The first one is that 

U_ = Ul (37) 
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Ml = Ml (38) 

i.e. + states and — states have equal mass eigenvalues but + states are mixed between 
them in a way just opposite to that of — states. As a consequence 

e_ = -e+ (39) 

Am^ = Am^ (40) 

E = (41) 
Inserting these in the expressions for the transition probabilities (p^-(p!9D we get 

P{^eL ^V^l) = (42) 

Pi'^eL^^tL) = sin2 2^+sin2^t (43) 

Pii^eL - i^Il) = (44) 

We see that, in this case, pure flavour oscillations and pure neutrino- antineutrino (Pon- 
tecorvo) oscillations are not allowed, while only flavour changing neutrino-antineutrino 
transitions are predicted. It is interesting to note that for the latter, the expression for 
the transition probability has the same form as (|33D . As we will remark below, this has 
implications on the interpretation of disappearance neutrino oscillation experiments. 

On the contrary to the ansatz just analyzed, we can further explore the possibility that 
neutrino mixing is given only by the Majorana mass term, while masses are generated 
directly by the Dirac term. We then consider the following mass matrices: 

Also in this case we find that relations (P7D,(|55|) hold, so that again we have eqs. 



(^) for the transition probabilities. The fact that both cases analyzed in this paragraph 
lead to the same phenomenological predictions is analogous to that encountered in the 
previous paragraph, where (|30D and (^) also gave identical results. We then observe a 
symmetry between Dirac and Majorana mass terms. 



2.3 The case of degenerate Dirac and Majorana mixing 

Another interesting ansatz is to suppose that neutrino mixing is generated by Dirac and 
Majorana mass terms with the same strength. This can be simply implemented by using 

Mn = ( It \ Mu=( ° \ (46) 



or 



Afo = f » \ Mm = ( < "J? ) (47) 
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In both cases we have one non-vanishing mixing angle and 2 (nearly) independent mass 
parameter; more precisely 

6^=0 (48) 
Am^ = cos20+Am^ (49) 



while S is given by 



2 sin^ 29^ (Am^'^ 



^ - (50) 
(mi_ + m2_) 

Note that in the limit of zero mixing angle we are left with only one mass parameter 
(S = 0). Instead for arbitrary mixing, in the present case, if the states z/i+, z/2+ are 
degenerate in mass (Am^ = 0) then a// the four states i'i±, i>2± are degenerate (Am^ = 
Ami = S = 0) so that no transition can occur. 

Inserting (|48|) , ( ^9|) in ([T7|)-(|T9D we now observe that all the transition probabilities are 
different from zero, and are given by 



1 Am^ 

Pil^eL ^iy,L) = P{yeL ^ V^l) = ^ siu^ 2^+ siu^ (51) 

9 ./S-2siAml\ „ . ./S + 2c+Ami\ 
Pi^eL^^L) = 4sin'(^ itj +4sm2(^ _± (52) 

- 44sm^(^^) (53) 
Interestingly, let us note that both Uei — > i^/iL and z/gL — > ^'^l transitions have the same 



probability, whose form is identical to ( pSf ) and (|43| ) apart a constant suppression factor for 
the oscillation amplitude of 1/4. For degenerate Dirac- Major ana mixing we then predict 
non vanishing flavour- conserving oscillations and equal probabilities for flavour- changing 
ones. 



3 Dirac- Major ana neutrino oscillations in matter 

In this section we generalize the MSW theory |^ to take into account the Dirac-Majorana 
nature of neutrinos. 

Let us consider neutrinos travelling in a medium with constant density, whose electron and 
neutron number densities are given by Ne and Nn, respectively. In the flavour basis, the 
evolution equations are simply given by (0) where the energy u is substituted hj uj — V,V 
being the effective potential experienced by a given neutrino state in the medium, namely 



26, 25 



K., = - K=„ = V2Gf (iV, - liV„) (54) 



eR 

G p 

V = Q for all the other sterile states (56) 



K,, = - K;^, = - ^ Ar„ (55) 
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Here Gp is the Fermi coupling constant, and for simplicity we consider only non mag- 
netized media (the generalization to these media is straightforward following ||2^). In 
compact form, the evolution equations are then given by 



where 



k + VL 
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In the ultrarelativistic limit, and for V^^^, V^^^ <^ fc, we can again separate the evolution 
of and ur states, and in particular we get 



HmTiL 

with the effective hamiltonian given by 
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(61) 



In terms of the physical mass and mixing parameters, subtracting from Hm terms propor- 
tional to the identity matrix which contribute with an irrelevant common phase factor to 
the wave functions, the effective hamiltonian in the flavour basis can be explicitly written 

as 

- a - D^r 
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where we have used the notations 
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Diagonalizing Hm in (|62D we then get the matter mass eigenstates from which the transi- 
tion probabilities in matter can be obtained. Before dealing with this point, we first want 
to discuss the occurrence of resonances in neutrino matter oscillations. 
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3.1 Resonance conditions 



Resonant enhancement of oscillations takes place when two unperturbed (mixing equals 
zero) energy levels cross between them The resonance conditions can then be (ap- 
proximately) obtained by equating the diagonal elements of Hm in (|6^). 
For resonant u^l —>■ i^^iL transitions we have the following condition 

— -i cos 26+ + -—- cos 26- = 2V2GFNe (66) 
while for Uei — > z^^l and U/^l i^Il 

cos 26+ + cos 26- = 2V2Gf (n, - ^ N^) (67) 



2k ^ 2k , y . 2 

Ami „ Am^ 



cos 26+ + — ; — cos 26- = V2GFNn (6^ 



2k ^ 2A; 



respectively. Instead, as already found in |]23|, the active-sterile flavour-conserving matter 



transitions VeL i^li ~^ i^^l have a maximum amplitude only if 

Vu, = (69) 



or 



Vu, = (70) 

respectively 5 



The relation (|66|) generalizes to Dirac-Majorana neutrinos the resonance condition for 

flavour oscillations obtained in [§] for pure Dirac or pure Majorana (ultrarelativistic) 
neutrinos: ^ 

cos26 = V2GFNe (71) 

An important novel feature emerging from (|66|) is that, on the contrary to what happens 
for pure Dirac or pure Majorana neutrinos (eq. (0)), the resonance condition is ruled by 
two squared masses differences so that, even if two or three of four mass eigenvalues are 
degenerate, the enhancement of oscillations can still take place (note that the resonance 
density in (|66|) is shifted with respect to that occurring for ( |7lD towards lower values). 
The same is valid for active-sterile (flavour-changing) transitions. Note, however, that the 
resonance conditions do not depend on the S parameter. 

For the particular case in which the Dirac and Majorana mass matrices are given 
by ( p6D or (p5D, only one transition {veL^ ^^.l) occur, and for this the resonance 
condition reduces to 

cos 26+ = V2Gf (n^ - (72) 



2k ' ^ V ^ 2 



^Note that, while eq. (|7^) is reaHzed only in vacuum, relation ( |69| ) can be satisfied also in a medium 
with Ne = Nn/'2.. Such a condition can be achieved in the first stages of the neutronization phase of a 
neutron star. 
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Instead, for the mass matrices taking the form in (^) or ( P7| ) all the transition discussed 
in the general case can be resonant, but now they are all ruled by only one squared masses 
difference. So, the resonance condition for — > t'/^L 

^^cos2^+ = V2GFNe (73) 

is phenomenologically equivalent to that for pure Dirac or pure Majorana neutrinos, while 
the one for z/g^. i^'^l again given by ([7^). 

After this qualitative discussion on the resonant enhancement of matter oscillations, 
we now proceed to find the expressions for the transition probabilities. As for the vacuum, 
this can be achieved also in the general case (described by the hamiltonian Hm in (|6^)) 
because the eigenvalue equation corresponds to a fourth degree algebraic equation whose 
solutions are known analytically. However, the resulting expressions for the probabilities 
are very involved and not very informative as (and much more than) in the vacuum case. 
So we will only consider the particular scenarios envisaged in the previous section and for 
these we will single out the explicit form of the transition probabilities. 



3.2 Pure Dirac and pure Majorana neutrinos 



Given (|30|) or (^) the problem is separable into the diagonalization of two 2x2 effective 
hamiltonians, one for the active left-handed neutrino states and the other for the sterile 
left-handed antineutrino states (which propagate freely as in vacuum) for the Dirac case 
or for the active right-handed antineutrino states for the Majorana case. In both cases 
we have 

Hml^L = ujul (74) 



with 



H„.^k + ^+{ '^^^ 1 (75) 




and so we recover the standard MSW theory 

No active-sterile (flavour-changing) transition occurs in this case, as in vacuum. 



3.3 The cases of Dirac mixing and Majorana masses and vice- 
versa 



Let us now assume that mass matrices are given by (pO]) or (^5]) and substitute relations 
(^)-(^) in the effective hamiltonian (^). Inspired from the results obtained for the 
vacuum case, it is useful to introduce the permutated flavour basis 



c 



\ 



V 



V. 



eL 



(76) 
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In this basis, the problem is again separable because the effective hamiltonian becomes 
block-diagonal. We then have 

(77) 
(78) 

with 

HS' ^t^^^^ilJ ) (80) 




where the mass matrix M. is 




^ = \ Jm Ji^ \ (81) 



Let us focus, for example, on (|7DD. It has the same form of the 2x2 effective e — MSW 
hamiltonian (^); the expressions for the transition probabilities in matter are then 

P{y.L ^ yuT,) = (82) 



\ oi,^^ ')Ul'^ oivi^ 



Pi^eL^^L) = sin2 2^!^sin2— (83) 



P{^eL ^Kl) = (84) 

(a; ~ t) where the effective mixing angle in matter is given by 

.in 2et = ^ " ^ (85) 

J(^ cos 29+ - ^Gf{N, - 1JV„)) + (^ sin 

while the effective oscillation length by 

27r 

Lm = , (86) 



cos 26+ - V2Gf{N, - iiV„)) + sin 2^^+) 



2 



From (pSD we then obtain again that the transitions u^l resonantly amplified if 

the resonance condition (|7^) is fulfilled. 

At this point we stress the fact that for Dirac-Majorana neutrinos described by the 
mass matrices in ( pO] ) or ( ^5]) the MSW theory [|] applies practically unmodified to 
^eL instead of u^l — ^ i^^l flavour transitions, and so the phenomenological impli- 

cations for disappearance experiments are the same for the two cases. In particular, the 
analysis performed in for the solar neutrino problem holds true for the present case, 
and from that the values of Ami and sin^ 29+ able to solve the puzzle can be extracted. 
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3.4 The case of degenerate Dirac-Majorana mixing 



Let us now turn on the last model considered in the previous section, with the mass 
matrices given by (^Bf ) or (^Tf) . From (^Sf) , (^Of) we then find that the proper effective 
hamiltonian to diagonalize is 



Hm = — 



1 

8k 



/ - 2 Ami cos 29^ 
— Ami 
E 

V - Am^ sin 26*+ 



— Am^ sin 29^ 
2 Ami cos 26'_^ 

— Ami sill 2^^_| 

S 



- Ami si'^ 26*+ 
2 Ami cos 26'^ 

- Ami sin 26'_i 



- Am^ sin 29+ \ 

E 

— Am^ sin 26*+ 
2 Am2 cos 29+ / 



Because of the presence of V^^^, we observe that the submatrices corresponding to the 
subsystem {ueL, ^^l) and {I'eL, ^"^l) ^'^^ equal (as instead happened for the vacuum 
case) so that UeL and Uei oscillations have different transition probabilities. 

The interactions with the medium remove this sort of degeneracy; this was manifest 
already in the expressions for the resonance conditions (|66D and (|67|). 
The exact eigenvalue problem for Hm in ( |87D involves again complicated solutions of 
a fourth-degree equation. Instead of giving the general expressions for the transition 
probabilities, it is more useful to discuss qualitatively the oscillations pattern which is very 
similar to the general one described by Hm in (^) since in both cases all the transition 
between the different flavour states can take place (the main differences between the two 
cases will be remarked in the following section). 

For a medium with varying density we can qualitatively analyze neutrino propagation in 
it with the help of a level crossing diagram, depicted (for the present case) in Fig. 1. Full 
lines represent the eigenvalues of Hm in (^^ plotted against the medium density, while 
the dashed lines correspond to the unperturbed (no mixing) energy eigenvalues (i.e. to 
the diagonal elements of Hm with 9+ = 0, since the non diagonal elements all vanish in 
the zero mixing limit). The crossing points of the unperturbed levels (approximatively) 



identify the resonance densities pm for the transitions z/eL — > z/^l(R1), u^l 



(R2) 



and VeL~^ ^liL (R-3). The relative positions of the three resonance points depend on the 
Yf. = Z/A ratio of the considered medium; as one can see from (|6^ 



PRl ~ K ^ 

1 -Y-' 



PR2 
PR3 



3K - 1 



-1 



(89) 

(90) 
= 0.5 



In Fig. 1 we have chosen Ye = 0.48; in the particularly common case in which Ye 
we would have that both R2 and R3 occur at the same density. 

Let us now focus on the evolution od Dirac-Majorana neutrinos propagating in a varying 
density medium, such as for example the Sun or another star like this. Neutrinos are 
produced deep in the star, at high densities, in flavour eigenstates (typically UeL ) and then 
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move out towards low density regions. As can be seen from Fig. 1, in fact, at high densities 
mixing effects can be neglected and the energy eigenvalues practically coincide with the 
unperturbed ones. The subsequent evolution, with decreasing density, depends whether 
the resonances are crossed adiabatically or not (i.e., qualitatively, if neutrinos travelling 
in the medium "see" a density which varies very slowly during their path or not). If the 
adiabaticity condition is fulfilled (at each resonance), neutrinos evolve according to 
the unbroken lines; otherwise at a given resonance there is a non vanishing probability for 
"jumping" from one level to another one p^. So, following Fig. 1, if for example we have 



a VeL — ^4 at high densities and the resonance R3 is crossed adiabatically, at intermediate 
densities we again encounter a but now 1/4 ~ i/^^ : we have had a transformation of an 
active v^l into a sterile z/^j^ . Instead, if R3 is crossed non adiabatically so that there is 
a non zero probability to jump onto the Vi state, at intermediate densities we encounter 
(with a certain probability) a z/i ~ and practically we have had no conversion. But 
moving towards lower density another resonance point (Rl) is present: if this is crossed 
adiabatically, then after that we have a z/i ~ u^l (active-active PeL^ i/^^ conversion) , 
otherwise we have a jump onto the V2 state which (after the resonance) is approximatively 
a VeL (no conversion). 

Confronting the present scenario with that of the usual MSW theory, we see that now 
efficient active-active flavour-changing conversions can be obtained only if the resonance 
for active-sterile transitions (R3 or R2) is crossed non adiabatically while that for active- 
active transition (Rl) is crossed adiabatically. But the important novel feature regarding 
Dirac-Majorana neutrinos is that exiting from the medium we do not have pure flavour 
states. In fact, as one can see from Fig. 1 and contrarily to what happens for the standard 
MSW theory, at very low density (approaching the vacuum) the energy eigenstates (solid 
lines) do not approach flavour eigenstates (dashed lines). This is a genuine feature of the 
Dirac-Majorana nature of neutrinos and it can be shown that (in the present case) this is 
due to the fact that in the limit of zero mixing the S parameter vanishes (see eq. (|50|) ) |. 
However, this is not a completely surprising feature, because we already know that in 
vacuum even for zero mixing the pure flavour eigenstates are not the physical eigenstates, 
which are instead given by the Majorana combinations n± in (|TB|). The level crossing 
diagram in Fig. 1 for zero density is just an expression of this physical fact. For zero 

mixing we have only two (doubly degenerate) energy eigenvalues ± while, switching 
on the mixing, four non degenerate energy eigenvalues appear, 

^(±2Am^^ + s) (91) 

^ (± 2Am^ cos 26+ - s) (92) 

and the physical eigenstates at the exit of the medium are just z/i, 1/2, ^3, = n± 
corresponding to the energy values given in (p^), (0). 



^In fact, keeping E fixed in the limit 0, at very low density the dashed lines approach the solid 

ones 
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4 Concluding remarks 



In this paper we have studied the propagation both in vacuum and in matter of Dirac- 
Majorana neutrinos and analyzed active-active (flavour-changing) as well as active-sterile 
transitions, which are, in general, both possible. 

For vacuum oscillations, in section 2 we have given the general expressions for the tran- 
sitions probabilities for PeL — ^ ^iiL , ^eL — ^ J^^l , ^ei — ^ t-'Il have then discussed some 
interesting limiting cases for Dirac (Md) and Majorana {Mm) mass matrices. For pure 
Dirac {Mm = 0) or pure Majorana {Mr, = 0) neutrinos obviously we recover the usual 
flavour oscillation formulae [0], while for both M^ and Mm non vanishing and diagonal 
the Pontecorvo formula for neutrino-antineutrino (active-sterile) oscillations is ob- 
tained p2 |. 

An interesting non trivial case is that with Md and Mm given by ( p6D or ( ^5] ) which imple- 
ment the idea that neutrino mixing is essentially ruled only by the Dirac mass term while 
the Majorana mass term is diagonal or vice- versa, respectively. In both cases, neither pure 
flavour oscillations nor Pontecorvo oscillations are predicted, but only flavour- changing 
active-sterile transitions, such as u^l — ^ i^^l ; possible. Remarkably, the transition 
probability for these is identical in form to that for flavour oscillations for pure Dirac 
or Majorana neutrinos, and this holds both in vacuum and in matter. For the latter, 
the resonance condition is only shifted by the neutral current contribution of u^l to the 
effective potential. So, for example, the solution to the solar neutrino problem in terms 
of active-sterile neutrino oscillations proposed in applies unmodified to the present 
scheme. 

Another interesting case, even if a bit more complicate, has been analysed for the mass 
matrices in (^6]) or (0), which implements the idea that neutrino mixing is given by 
the Dirac and Majorana mass terms with the same strength. In this case, ^iiL , 

^eL 5 ^eL ^eL transitions are all possible and, in vacuum, the first two have 

the same transition probability, which is also identical in form to that obtained in the 
previous case, except for a constant suppression factor in the amplitude of oscillations 
of 1/4. Also in matter the pattern of neutrino transitions present in the general case is 
(qualitatively) reproduced in this peculiar scheme. In particular, all the flavour changing 
oscillations can be resonantly amplified while Pontecorvo — > i^ei matter oscillations 
have maximum amplitude only for a given electron to neutron number ratio (see eq. (|69|) 
and the related footnote); the resonance conditions were discussed in section 3.1 . 
Given the multiresonance structure of the oscillations pattern, it is then interesting to 
follow the evolution of a , for example, in a varying density medium such as the Sun; 
this has been done in section 3.4 with the help of a level crossing diagram reported in Fig. 
1. Several scenarios are possible according to the adiabaticity properties of level crossing 
near the resonance points. In particular, starting from a pure u^l beam at high density, 
to have a consistent conversion into at low density the resonance for u^l — > t^^l to 
be crossed non adiabatically, while the passage through the one for u^l — > i^^il has to be 
adiabatic. However, we have also shown that at very low density, and then in the vacuum, 
it is more appropriate to deal with the Majorana combinations n± in (pTSf ) than with the 
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pure flavour states v^l , ^eh i ^p-l ; ^^^l • This is strictly related to the Dirac-Majorana 
nature of neutrinos, which chooses Majorana eigenstates instead of pure flavour states 
as starting points. In this respect, we have to deal with "generic" flavour-changing or 
flavour-conserving transitions of Dirac-Majorana neutrinos without looking at the partic- 
ular active neutrino or sterile antineutrino state. It is through the weak interactions, with 
which neutrinos are produced and detected, that a particular (active or sterile) compo- 
nent of the Majorana eigenstates is chosen. 

The results obtained for the case of degenerate Dirac-Majorana mixing are qualitatively 
valid also in the general case in which all the entries of the Dirac and Majorana mass 
matrices are non zero and different between them. The main difference between the two 
cases is that in the general framework there are 3 mass parameters and two mixing angles 
ruling the evolution, while for the particular case studied in sections 2.3 and 3.4 there 
are only 2 mass parameters and 1 mixing angle (these parameters being not completely 
independent, because of relation (|5DD). The presence of more degrees of freedom in the 
general case allows to consider some peculiar situations which are not possible otherwise. 
The most remarkable one is that in the general case the proportionality of the S parame- 
ter to sin^ 26'+ (see eq. (|3)) is lost, so that the structure of the level crossing diagram at 
very low density can be altered. The eigenvalues of Hm in (|62D for zero density (vacuum) 
are given by 

^ (±2Am^ + S) (93) 

^ (±2Am2_ - S) (94) 

so that one can manipulate the mass parameters to modify the low density region of the 
level crossing diagram without grossly altering the region where the resonance points are 
present. In any case, there can be present no substantial modifications of the conclusions 
reached above. 

The oscillations of Dirac-Majorana neutrinos here studied with their peculiar features 
can be efficiently tested in astrophysics, in particular detecting solar or supernova neu- 
trinos, and can have even profound implications in cosmology for the nucleosynthesis of 
light elements in the Universe. 
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Figure 1: Level crossing diagram for 1 MeV momentum Dirac-Majorana neutrinos 
described by the mass matrices in (H) or (^T]) . The eigenvalues of Hm in (|87|) are plotted 
versus the density of a medium with Ye = Z/A = 0.48. Neutrinos parameters are fixed 
as follows: Ami = ^O'^eV^, E = 2.5 x 10"^ e1/^ sin 26+ = 0.1. Dashed lines refer to 
the zero mixing limit {6+ = 0) of the energy eigenvalues. 
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